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Abstract 

We study the spectrum of the QCD Dirac operator by means of the valence quark 
mass dependence of the chiral condensate in partially quenched Chiral Perturbation The- 
ory (pqChPT) in the supersymmetric formulation of Bernard and Golterman. We consider 
valence quark masses both in the ergodic domain (m^, <^ Ec) and the diffusive domain 
(m„ ^ Ec). These domains are separated by a mass scale Ec ~ F'^/'EqL^ (with F the 
pion decay constant, Sq the chiral condensate and L the size of the box). In the ergodic 
domain the effective super-Lagrangian reproduces the microscopic spectral density of chi- 
ral Random Matrix Theory (chRMT). We obtain a natural explanation of Damgaard's 
relation between the spectral density and the finite volume partition function with two 
additional flavors. We argue that in the ergodic domain the natural measure for the su- 
perunitary integration in the pqChPT partition function is noncompact. We find that 
the tail of the two-point spectral correlation function derived from pqChPT agrees with 
the chRMT result in the ergodic domain. In the diffusive domain we extend the results 
for the slope of the Dirac spectrum first obtained by Smilga and Stern. We find that the 
spectral density diverges logarithmically for nonzero topological susceptibility. We study 
the transition between the ergodic and the diffusive domain and identify a range where 
chRMT and pqChPT coincide. 



PACS: 11.30.Rd, 12.39.Fe, 12.38.Lg, 71.30.+h// /noindent Keywords: QCD Dirac opera- 
tor; Chiral random matrix theory; Partially quenched chiral perturbation theory; Thouless 
energy; Microscopic spectral density; Valence quark mass dependence 



1. Introduction 



Spontaneous breaking of chiral symmetry is believed to be an important property of 
the strong interaction. In the hmit of massless quarks the QCD Lagrangian is invariant 
under UL{Nf) x UnlNf) transformations, but the vacuum state is not. The analysis of 
hadron spectra and numerical simulations ^ ^ on the lattice strongly support this 
assertion. The order parameter for spontaneous breaking of chiral symmetry is the chiral 
condensate. In the chiral limit it is related to the spectral density of the QCD Dirac 
operator by means of the Banks-Casher formula 0], 

S„ = l.m^. (1) 

where V is the volume of space-time and the spectral density is defined by 

p(A)=^5(A-A,). (2) 

k 

Here, the are the eigenvalues of the Euclidean QCD Dirac operator. They will be called 
virtualities in this paper. In (|l]) the thermodynamic limit is taken before the chiral limit. 
For broken chiral symmetry, the eigenvalues near zero are spaced as l/p(0) = it/J]qV. In 
order to study the behavior of the Dirac spectrum in the approach to the thermodynamic 
limit it is therefore natural to introduce the microscopic variable |^ 

u = VXEq, (3) 

and the microscopic spectral density defined by 

Psiu) = hm T^ipiv^))- (4) 

In general, the average (■ ■ ■) in this equation is over the probability distribution of the 
eigenvalues. It was conjectured in that the microscopic spectral density is universal and 
is given by chiral Random Matrix Theories (chRMT) with the global symmetries of the 
QCD partition function. This conjecture has been verified both by instanton liquid and 
by lattice QCD |]^, ||, ^ |T0| simulations. It was also found from lattice QCD simulations 
that eigenvalue correlations on the microscopic scale, both near zero virtuality and in the 
bulk of the spectrum, are given by Random Matrix Theory [^, Q . The foundation 

in 



for these observations was laid by universality studies [[T^, |T5|, [T^, |T9|, ^ 

which it was shown that the microscopic spectral density and other spectral correlations on 
the scale of individual level spacings are not affected by macroscopic deformations of the 
random matrix ensemble. The interpretation is that quarks undergo a chaotic /diffusive 
motion in four Euclidean dimensions and one artificial time dimension. Such picture is 



also suggested by the success of microcanonical lattice QCD simulations p2[- According 



2 



to the Bohigas conjecture the correlations of the eigenvalues on the microscopic scale 
are then given by Random Matrix Theory. In Ref. [24| chiral symmetry breaking was 
also related to transport properties of quarks in 4+1 dimensions. 

Of course, the dynamics of QCD, which is not included in chRMT, should enter in 
the Dirac spectrum. Therefore, we expect a scale above which the microscopic spectral 
density and the spectral correlations are not given by chRMT anymore. The existence of 



such scale is well known in the context of mesoscopic physics (seep5|, for reviews). 

It is known as the Thouless energy, which is the inverse diffusion time of an electron 
through the sample 



hD 



(5) 



where D is the diffusion constant and L is the linear dimension of the box. In the case 
of QCD, the diffusion process is that of a quark propagating through the Yang-Mills 
fields in a 4+1 dimensional space time. A second energy scale that enters in mesoscopic 
physics is ^/xg, where is the elastic diffusion time. Based one these two scales three 
different domains for the energy difference, that enters in the two-point correlation 
function, can be distinguished: the ergodic domain, the diffusive domain or Altshuler- 
Shklovskii domain and the ballistic domain. For energy differences 5E <C E^ eigenvalue 
correlations are given by RMT. Since for time scales larger than the diffusion time, an 
initially localized wave packet explores the complete phase space, this regime is known 
as the ergodic regime. In the diffusive regime, defined by E^. <^ 5E <^ Ti/Te only part 
of the phase space is explored by an initially localized wave packet, resulting in the 
disappearance of eigenvalue correlations. In this paper we don't consider the ballistic 
regime with 5E ^ Ti/re- For an interpretation of the Thouless energy in terms of the 
spreading width we refer to pOl As was shown in ||3l|, |3^, the spectral two-point 



function can be related to the semiclassical return probability which provides a simple 
intuitive picture of its asymptotic behavior. For other recent studies on this topic we refer 
to 1^ |3^, ^ |33, ^ . What has emerged from these studies is that there is a close 
relation between eigenvalue correlations and localization properties of the wave functions. 
The interpretation of spontaneous chiral symmetry breaking as a delocalization transition 
was made earlier in |^0|. By analogy with the Kubo formula, Sq plays the role of the 
conductivity 

In QCD Dirac spectra the Thouless energy was identified as the scale where the inverse 
pion mass is of the order of the linear size of the box |^, Using standard relations 
based on the approximate chiral symmetry of the QCD partition function it has been 



suggested that |4l|, ||, |2 



Er. 



(6) 
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where F is the pion decay constant. The existence of such a scale has been observed in 



instanton hquid |41, H3] and lattice QCD M, R[ simulations. One of the main goals of this 



paper is to identify this scale by means of partially quenched Chiral Perturbation Theory 
In this paper we will discuss four different partition functions: the Euclidean QCD 



45 



partition function, the effective finite volume partition function [46, the chiral random 



matrix theory partition function [E[ and the partition function for partially quenched 



Chiral Perturbation Theory |i5]. Below, we briefly define the first three partition functions 
whereas the last one is discussed extensively in the remainder of this paper. At this point 
we only mention that pqChPT is an effective mesonic theory which is generally believed 
to describe the low-energy limit of QCD. In addition to the usual Nf sea quarks, this 
theory has k (in this paper k = 1) valence quarks and an equal number of bosonic ghost 
quarks with the same mass rriv Without a source term that breaks the supersymmetry 
the valence sector drops out of the partition function. In other words, the valence quark 
mass only appears in the operators. 

The Euclidean QCD partition function with fundamental quarks is given by 

Zqcd{M) = JdA det(^ + M) e-^^^, (7) 

where is the quark mass matrix, I/) is the Euclidean Dirac operator and Sym is the 
Yang-Mills action. The integral is over all gauge field configurations A. The eigenvalues 
of the Dirac operator are defined by 

Iptpk = i^k'ipk- (8) 

Notice that for 7^ the eigenfunctions occur in pairs iIju, 'J^i'k with opposite eigenvalues. 
The quark mass dependence resides in the determinant. 

For an Euclidean space-time volume with linear dimension in the range 

l/A<L<l/m^, (9) 

(A is a typical hadronic scale and m-„ is the pion mass) the mass dependence of the QCD 
partition function is given by the effective finite volume partition function. For Nf flavors, 
vacuum angle 9 and mass matrix JH it is defined by [|^, ^ 



Zf(M,9)r^f rf[/e^^o^^^^^^"'"^ (10) 
JueG/H 

The integral is over the Goldstone manifold associated with chiral symmetry breaking 
from G to H. For three or more colors with fundamental fermions G/H = SUlINj) x 
SU{i{Nf)/ SUA{Nf). The effective finite volume partition function in the sector of topo- 
logical charge u follows by Fourier inversion with respect to vacuum angle 6. The partition 
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function for z/ = is thus given by (|ToD with the integration over SU{Nf) replaced by 
an integral over U{Nf). For an elaborate discussion of this partition function we refer to 



47 



The chiral random matrix partition function with the global symmetries of the QCD 
partition function as input is defined by ^ 



Nf 

Z^^_,(mi,---,m^^)= |wndet(P + m;)e-^^^^^^, (11) 



where 



/=i 



and W is a. n X m matrix with u = \n — ml and N = n + m. As is the case in QCD, 
we assume that the equivalent of the topological charge u does not exceed V^, so that, 
to a good approximation, n = N/2. Then the parameter E can be identified as the 
chiral condensate and as the dimensionless volume of space time (Our units are defined 
such that the density of the modes N/V = 1). The matrix elements of W are either real 
{(3 = 1, chiral Gaussian Orthogonal Ensemble (chGOE)), complex (/? = 2, chiral Gaussian 
Unitary Ensemble (chGUE)), or quaternion real [jS = 4, chiral Gaussian Symplectic 
Ensemble (chGSE)). For QCD with three or more colors and quarks in the fundamental 
representation the matrix elements of the Dirac operator are complex and we have (3 = 2. 
It can be shown that in the domain (P) the random matrix partition function can be 
mapped onto the effective finite volume partition function p . For more discussion of the 
chRMT partition function we refer to ||49|| . 

In this paper we study the Dirac spectrum by means of the valence quark mass de- 



pendence of the chiral condensate [50, 51 



where (■ ■ ■) denotes an average with respect to the distribution of the eigenvalues. The 
spectral density follows from the discontinuity across the imaginary axis. 



DiscL^^,,S(m,) = limS(zA + e) - T.{i\ - e) = 2nJ2{6{X + A,)) = 27r(p(A)). 



(14) 



The two point spectral correlation function is given by 

— Disci 

4^2 lm„=iA,m„,=a' ^ \ + iXi + m^, 



(P(^)P(^')) = DiscL„^,,^^,^,, E ( a.L,.a,LJ • ^15) 
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On the other hand, the valence quark mass dependence can be calculated from the par- 
tially quenched partition function. By differentiating with respect to the valence quark 
masses we obtain in the parametrization discussed in section 2, 



1 /■ /■ 

(p(A)p(A'))^ = ^ Disc|„^^.;^ ,„^,^.;^, J d^x J 



(16) 



where the superscript C refers to the connected correlation function which is defined as 
(p(Ai)p(A2))'^ = (p(Ai)p(A2)) - (p(Ai))(p(A2)). The notation {■■■)pq denotes averaging 
with respect to the partially quenched partition function. The meson fields with quark 
content corresponding to a mass M = {m^ + myi)T,Q/F^ are denoted by ^^'^ The expec- 
tation value of the product of the meson fields can be expressed in terms of the meson 
propagator in momentum space. Taking into account the combinatorial factors and per- 
forming the integrals over x and y and one momentum integral this results in 

1 1 

(p(A)p(AO)^ = — DiscL„=a,™„,=iA' E ^q2^M^y ^^^^ 

We have expressed the two-point correlation function as a sum over momentum modes. 
Such relations are well known from the theory of disordered systems (see for example 



33| , |31| , p^). The zero momentum modes reproduce the asymptotic limit of the random 
matrix result for the two-point correlation function. This can be seen by evaluating the 
discontinuities after substitution of the Gell-Mann-Oakes-Renner relation for the meson 
mass. For |A — A'| ^ A but much larger than the average level spacing we obtain 

(P(A)P(A')>- ~ -i^x^. (18) 

which is the well-known random matrix result. The full correlation function can only 
be obtained by means of an exact integration over the saddle-point manifold |5^. In 



this paper we will not perform this task for the two-point function but instead focus on 
the one-point function. Because of the Ua{^) symmetry, the microscopic spectral density 
already shows a nontrivial structure for eigenvalues close to zero as has been observed in 
lattice QCD simulations 0, ||. The discussion of the two-point function will appear in a 
future publication p^ . 

The valence quark mass dependence of the chiral condensate was calculated in par- 
tially quenched chiral perturbation theory by Golterman and Leung [E^] by means of 



a supersymmetric method introduced by Bernard and Golterman [0, We use this 
result to calculate the spectral density in the perturbative domain. Among others we 
reproduce the slope of the spectral density first derived by Smilga and Stern . When, 
for virtualities below the Thouless energy, the zero modes are treated exactly, we repro- 
duce the random matrix result for the microscopic spectral density. In fact, we show that 
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the partially quenched supersymmetric partition function can be written as an ordinary 
effective finite volume partition function with two more flavors which are the equivalent of 
one valence quark flavor and its supersymmetric partner. In this way we obtain a natural 
explanation of a result first found by Damgaard \^% . 

The organization of this paper is as follows. Partially quenched chiral perturbation 
theory is introduced in section 2. In section 3 we discuss its different domains of validity. 
The choice of the superunitary measure is motivated in section 4. In sections 5 and 6 the 
partially quenched partition function is studied in the domain where it is dominated by 
the zero modes. The valence quark mass dependence in the quenched limit is calculated in 
section 5, and, in section 6, we show that the microscopic spectral density that follows from 
the valence quark mass dependence of the chiral condensate coincides with the result from 
chRMT. In section 7 we calculate the valence quark mass dependence of the condensate by 
means on pqChPT. The transition from chRMT to pqChPT is discussed in section 8. In 
section 9 we evaluate the spectral density in the diffusive domain. Concluding remarks are 
made in section 10. An explicit calculation of the spectral density for Nj = 1 is presented 
in Appendix A and a short derivation of a Berezinian that occurs in the diagonalization 
of superunitary transformations is given in Appendix B. 



2. Partially Quenched Chiral Perturbation Theory 

Partially Quenched Chiral Perturbation Theory (pqChPT) is an effective theory de- 
scribing partially quenched QCD at low energies |^5|, This is a theory with both sea 
quarks, appearing in the fermion determinant, and valence quarks which only appear in 
the operators. Partially quenched chiral perturbation theory interpolates between stan- 
dard (unquenched) Chiral Perturbation Theory (ChPT) |Q and completely quenched 
Chiral Perturbation Theory (qChPT) [6^, |61| . In this section we give a brief overview of 
partially quenched chiral perturbation theory in the formulation of Bernard and Goiter- 
man. For more details we refer to the original literature 

We consider a QCD-like theory with one valence quark and its supersymmetric partner, 
both with mass m^, and Nf unquenched quarks (sea quarks) with mass m^. Then the 
fermion determinant corresponding to the valence quark is cancelled by the contribution 
from a bosonic ghost quark with the same mass |Q. The full chiral symmetry of the 
theory is thus the direct product of the supergroups (also called graded groups) 

G = UL{Nf + l\l)(^UR{Nf + l\l), 

where one of the U{1) groups is broken by the anomaly. The partially quenched QCD 
Lagrangian at a given order in the momentum can be approximated by an effective La- 
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grangian for the field U G SU{Nf + 1|1) U{1) which transforms linearly under the chiral 
group. (Here, denotes a semi-direct product.) In terms of the Goldstone fields U is 
parametrized as 

U = exp(iv^$/F), (19) 
where $ is a {Nf + 2) x [Nf + 2) super- Hermit ian matrix field: 

and F is the pion decay constant. Here, is the {Nf + 1) x [Nf + 1) matrix containing the 
ordinary mesons made of quarks and antiquarks, (p is the meson with two ghost quarks, 
and which is a vector of length Nf + 1, represents the fermionic mesons consisting 
of a ghost-quark and an ordinary anti-quark. We choose a diagonal mass matrix with 
degenerate sea quark masses resulting in the [Nf -|- 2) x [Nf + 2) quark-mass matrix 

M = dmg[my,ms, . . . ,ms,m^). (21) 

The corresponding pion masses will be denoted by 

M^. = [rrii + mj)So/F2 for z,j = s, v. (22) 
Expanding the Euclidean effective Lagrangian in terms of the momenta and the quark 



masses we obtain to lowest order ||4^, |5^ 

Aff = ^ StT[d^U d,U^) - ^ Str[MU + MU^) + ^ $ ^ + 5 d,^o d,^o, (23) 
4 2 

where we have added a mass term and a kinetic term for the super-r/' flavor-singlet field 
$0 = Str($) with parameters ttiq and a (for the definition of the supertrace, Str, (also 
called graded trace) we refer to [^). In the next section we will show that the term 
containing a is sub-leading in the chiral expansion. The quark condensate in the chiral 
limit is denoted by Sq. The partition function corresponding to (^) is defined by 

Zpg = J rff/e-/'^'"^^ff (24) 

and will be studied for different ranges of the quark masses. 

A physical interpretation of the zero momentum component of the flavor singlet mass 
term can be obtained by linearizing the square according to 

2 1 » 2 



^-v^% ^ / duSdet^U) e"^, (25) 
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where 



(^^> ^ ^ (26) 

is the famihar topological susceptibility. For the definition of the superdeterminant, Sdet, 
(also called graded determinant) we refer to This result allows us to interpret ml 



as a topological contribution to the singlet mass. The partition function with zero total 
topological charge is obtained by multiplying the Goldstone fields with the Ua^^) phase 
and extending the integration manifold to include this contribution (see |^). In the 
zero momentum sector, the partition function for u = is thus given by (|2^ ) and the 
Lagrangian (^) without the topological mass term mQ$g/6 but with the integration 
domain extended to U{Nf + 1|1). This partition function will be discussed in sections 
5 and 6. Notice that the microscopic spectral density depends on the global topological 
charge. On the other hand, the spectral density at a macroscopic distance from zero is 
not sensitive to the global topological charge. 



3. Domains of Validity 



For standard chiral perturbation theory the different domains of its validity were dis- 
cussed in detail by Gasser and Leutwyler . Since the discussion for partially quenched 
chiral perturbation theory is quite analogous we refer to this paper for additional argu- 
ments. 

With the basic scale given by e = 1/L <C ^qcd, we consider five different mass 
scales: the smallest eigenvalue of the Dirac operator Amin ~ e^, the valence quark mass 

~ e", the sea quark mass ~ e^, the pion momenta ~ 1/L ~ e and the 1]' mass 
■rrijji ~ Aqcd ~ 1- 

The integral over U in can be split into an integral over the zero momentum 
modes, Uq, and the nonzero momentum modes, ^(x), i.e. 



U = Uoe 



(27) 



The two different types of modes are mixed by the mass terms and by the singlet in- 
teraction. Keeping only terms up to second order in ^, our Lagrangian can be written 
as 

So 



d xC, 



eflf 



Str 



+ 



J d'x (^^d,^, d,ia + !^ $2 ^ ^ d,^^ . (28) 
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The nonzero mode fields are 0(e) whereas the zero mode fields are 0(e*^^~°^/^). The term 
that mixes both types of modes is of order independent of the value of a (for 1 < a < 4). 
This implies that the mixing term is always subleading. Since the fluctuations of the 
singlet fleld are of O(e^), the kinetic term of the singlet fleld is subleading as well. 

Since the mass term of the nonzero modes is of order e"~^ the integral over the nonzero 
modes does not contribute to the mass dependence of the partition function for a > 2. 
Therefore, for a < 2 (or m < the mass dependence is completely determined 

by the zero modes. The crossover-point can be determined more accurately by considering 
the mass dependence of the propagator of the nonzero modes. We will return to this point 
in section 8. 

In the chiral limit, —* 0, the integrals over the zero momentum modes become 
perturbative for ^ ^min, whereas for -C the mass dependence of the 

partition function is completely determined by the zero momentum modes. In the domain 

Amin < < 1/Aqcd^^ (29) 

the valence quark mass dependence of the chiral condensate can thus be determined both 
from the integral over the zero momentum modes and from a perturbative calculation. In 
section 8 it will be shown that this consistency condition is satisfled. 

Actually, because of the mixing by the singlet mass term, the above discussion has to 
be reflned by considering modes that are eigenvectors of the mass matrix. This is essential 
for the quenched limit. In this case the relevant mesonic mass matrix in a quark basis is 
given by |]60[ 



3 

3 -"'■'dd 



(30) 



To lowest order in M^^ the eigenvalues of this matrix are given by 2ml/3 + M^^/ (2mQ/3) 
and — M^„/(2mg/3). The zero momentum mode corresponding to the second eigenvalue 
only becomes perturbative for ^ e^, where also the nonzero modes contribute to the 
mass dependence. We conclude that in the quenched case there is no common domain of 
validity. 

Let us contrast this with the mesonic mass matrix for Nf = 1 which is given by 



-2 " 

^2 3 



V 



^~3~ 3~ ~3~ 



2 


\ 






3 




vv 


) 



(31) 



In the chiral limit (M^s — ^ 0) the eigenvalues up to order M^^ are given by Xi = ml + 
2M^J3, As = -M^JVS - M^J3 and A3 = M^jV^ - M^J3. In this case, the zero 
momentum modes corresponding to the smallest eigenvalue become perturbative for rriv ^ 
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Xmin and a common domain of validity can be identified. For Nf massless sea quark flavors, 
Nf — 1 eigenvalues are zero in the diagonal basis. However, mesons with two sea quarks 
do not enter in the one-loop calculation of the valence quark mass dependence of the 
chiral condensate. The smallest non-zero eigenvalue in the valence sector is again of order 
M^^. Therefore, as for one flavor, a common domain of validity can be found. Of course, 
if also the sea quark masses are in the range (j2^) a common domain of validity for the 
dependence of the chiral condensate on the sea quark masses can be found as well. 



4. Choice of Measure 



For the group integration in (123), we have two different possible choices for the supe- 
runitary integration measure: a compact one which was used in ||60| , ^ , and a noncompact 

I in the context of chiral Random Matrix Theory. 



one which was introduced in 53 



It was argued in JT^ that a noncompact measure reproduces the invariance of the bosonic 
degrees of freedom under noncompact Ua{^) transformations and therefore should be the 
appropriate measure for the quenched theory. We stress that for perturbative calculations, 
both measures lead to the correct result. 



as m 



60, 45 



For the purely quenched theory a decomposition of the compact measure in terms of 
the eigenvalues of the superunitary matrix U can be written down easily. The compact 
measure for the group can thus be written as 



DU 



e'" - e 



id>\2 



DU' 



in the parametrization 



U = U'AU'-\ 



(32) 



(33) 



with A a diagonal matrix with eigenvalues exp(z6') and exp(i0). This decomposition 
can be easily generalized to an arbitrary superunitary group. The generalization of the 
Berezinian appearing in (|3^ is given in Appendix B. 

The measure for the noncompact version of t/(l|l) is obtained by analytic continuation 
to imaginary angle 6 —is with s G (— oo; oo). This results in the measure 

dsd(f) 



DU 



DU' 



(34) 



(e^ - e''t'){e-'' - e-*<^, 

which coincides with the measure that was obtained in after a lengthy calculation. 

We will show in section 5 that the quenched partition function with this noncompact 
measure reproduces the quenched valence quark mass dependence of chiral random matrix 
theory. 
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The generalization of this measure to an arbitrary superunitary group appears to be 
more involved than this naive analytic continuation. However, it can be shown that 
the integration contours contributing to the discontinuity of the valence quark mass 
dependence of the chiral condensate across the imaginary axis can be deformed such that 



only an integral over a compact domain remains |5^, ^ 0. In terms of the partition 
function this corresponds to replacing the noncompact measure by a compact one. This 
will be the starting point of the calculation of the microscopic spectral density discussed in 
section 6. The calculation of the valence quark mass dependence of the chiral condensate 
for arbitrary Nj using a noncompact measure will be discussed elsewhere ||64 . 



5. Valence Quark Mass Dependence of the Chiral Condensate in the quenched 
limit 



In this section we will show that in the quenched case the noncompact measure (|3^) 
leads to the correct result for the valence quark mass dependence of the chiral condensate. 
As argued at the end of section 2, the partially quenched partition function in the sector 
of zero topological charge is obtained by extending the integration over U in (p^) to 
U{Nf + 1|1) and ignoring the singlet channel mass term in (]23|). In the quenched limit the 
integration is thus over ?7(1|1). With the measure ( p^ the sector of the zero momentum 
modes of the partition function (E^ in the quenched limit can be written as 



Z{J) 



dUdOds 



exp 




cos^^ 







cosh s 



u- 



The normalization of this partition function, 2'(0), follows from Wegner's theorem for 
superinvariant functions The integral over U can be performed by means of the 

|66|, |67l|. After differentiation 



supersymmetric version of the Itzykson-Zuber integral 
with respect to J the result for the valence quark mass dependence of the chiral condensate 
is 

cos 9 — cosh s 



— I dOds- 



7T 



^^(cos 6— cosh s) 



gs — e*^)(e * — e 



(36) 



where we have introduced the notation x = rnVTiQ. This integral can be evaluated most 
conveniently after differentiation with respect to x 



So 

4:71 



dOdsil 



)(1 



id+s\ xicos 6— cosh s) 



(37) 
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where the integrand has been simphfied by means of the identity 

cos 9 — cosh s 1 



e 



ie _ gs 2 



l-e-^e-'"). (3^ 



The integrals appearing in (^) factorize and can be expressed in terms of the modified 
Bessel functions Kq and Iq, 

d.Eix) = ^o[Io{x)Ko{x) - h{x)K^{x)]. (39) 

The r.h.s. of this equation can be rewritten as dx[xIo{x)Ko{x) + xIi{x)Ki{x)]. This 
equation can be integrated trivially with an integration constant that follows from the 
asymptotic limit S(x) ^ So for x — > oo. The final result is 

S(x) = Sox[/o(x)iro(x) + hix)Ki{x)], (40) 

which agrees with the valence quark mass dependence obtained from chiral Random Ma- 
trix Theory |5l|. In its domain of validity, this result reproduces the valence quark mass 
dependence found in lattice QCD [O, E^. 



6. Microscopic Spectral Density in the Ergodic Domain 

In this section we consider the zero momentum component of the partially quenched 
partition function for zero total topological charge. In it was shown that the 

noncompact integral contributing to the imaginary part of the chRMT partition function 
can be rewritten as an integral over [0,27r]. The spectral density can then be derived 
from a compactified partition function. By a direct translation of these arguments to the 
partition function of pqChPT, we obtain for the spectral density 



1 1 (97° 

(p(A)) = — Im ^« 



2tt Z0(J = 0) dJ 



(41) 

=a,j=o 



where [J) is the zero momentum pqChPT partition function with the com|)aci measure 



defined by 

Z'(M) = f rft/e^s*^(^^+*"^'). (42) 



Here, the integral is over ?7(A^/ + 1|1) corresponding to the sector of zero topological charge. 
The mass-matrix M is diagonal with matrix elements (m^ + J, m^, ■ ■ ■ , mg, — J). 

'4, for ^/ 



In a diagonal representation of the of superunitary matrices f/, for Nr flavors can 



be written as 

Z^M) = J n^^4c^t/|S(A)re^s*^*^^(^+^*)^"\ (43) 
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Here, |-B(A)p is the Berezinian for the transformation to eigenvalues and eigenvectors as 
new integration variables (see Appendix B). The quantity B{A) is given by 

= n,.(Ai-M) ■ 

where the bosonic and fermionic eigenvalues of U are denoted by and A{, respectively. 
In our case we have Nf + 1 fermionic eigenvalues and one bosonic eigenvalue. In the 
quenched case we use the supersymmetric normalization with Zp^{J = 0) = 1 whereas 
in the partially quenched case Z^^^J = 0) = (defined in eq. (p!OD). Here, and below 
the constants multiplying the partition function are absorbed in the integration measure. 
The integration over U G U{Nf + 1|1) can be performed by means of a supersymmetric 
version of the Itzykson-Zuber integral p^, |60, p^. For J ^ this results in 



P9V ; / 11 " B{cosek)B{MkVJ:o) o<k,i<Nf ^ ^ 



k 



For J = 0, the Grassmannian components of the integration paths lead to additional con- 
tributions. These so called Efetov-Wegner |69| terms only enter in the normalization 



of the partition function (see [[TOl , p^ , pO| , [7T[| for more discussion of this point). 

We evaluate this integral in the chiral limit, 0, for a fixed value of the valence 

quark mass. In this limit the ratio of the determinant and B[MkVT,o) is finite. In 
Appendix A, we evaluate the integral for Nf = 1 and indicate how it can be calculated 
for arbitrary Nf. In the remainder of this section we show that this generating function 
can be expressed as an ordinary effective finite volume partition function (see (|10|)) with 
Nf + 2 flavors. 

For arbitrary Nf the partition function ( |45| ) can be written as 



with fiQ = {niy + J)VT,Q, ^Nf+i = {^v — J)y^o and /i^ = m^VSo, k = 1,- ■ ■ Nf. For 
equal sea quark masses the ratio of the last factor and I\k<i{f''k — fi'i) has to be understood 
as a limit. The normalization factor that enters in the partition function is denoted by 
J\f. The function uj{9, 0) defined by 



uj{e,<j>) = ^n 4, (47) 

COS — COS 



has the important property that 



, ^ , , 4 , cos 6 — cos (h , , 
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This allows us to bring all factors uj{9k-,0i) into the numerator. Slightly rearranging the 
different factors and changing the sign of fiNf+i we can write the partition function in the 
more symmetric form 



cos t 



0<k<Nf k 0<k<l<Nf+l yr'k H'l) ^^/"T^ 

(49) 

where the diagonal mass matrix M is defined by M = {m^ + J, m^, ■ ■ ■ , m^, —m^ + J). In 
this integral we recognize the factors that occur in the ordinary Itzykson-Zuber integral. 
It can thus be rewritten as 

Z'„{M) = AT' n ifil - /.^^^J / dUl[d9k\A{e^'^)\'e'^^'^^^-^'^^'''\ (50) 

0<k<Nf •' k 

where the ?7- integral is over U {Nf + 2). Up to a constant factor this is exactly the effective 
finite volume partition function for Nf + 2 flavors. In the microscopic limit with quark 
masses of order ~ l/\^So, this partition function is equivalent to the random matrix 
partition function with Nf + 2 flavors. For masses in the microscopic domain we thus 
have 

Zl^{M)=U" n {^4-^^l,^^)Z^RJ!,i^\m, + J,m,,...,mM,,-m, + J), (51) 

0<k<Nf 

where Z^^j^'!^^'^ is the chRMT partition function as defined in (|ll]) for Nf + 2 flavors, zero 
topological charge and (3 = 2. Since the RMT partition function is even in all masses (|5T|) 
can be rewritten as 

Zlq{M) = U" n (^fc - l^lj;+i)ZmiT^\T^v + J, mi, ...,mNf,m^ - J). (52) 

0<k<Nf 

If we notice that /ig — fJ^'j^^.^i = 4mi,SoFJ + 0( J^), the spectral density is easily obtained 
from (HII) by differentiating with respect to J and setting rriv = iX- After taking the 
imaginary part we obtain for the microscopic limit of the spectral density defined in (|^), 

y{Nf+2),. . X 

Ps{u) = -\u\ [[ ifik + u) -jv— , (53) 



l<fc<Af 



where u = WHq. The denominator arises from the normalization Zp^{J = 0) of the 
partially quenched partition function. The overall normalization is such that p(A) ^ Sq/vt 
for A — > oo. This result coincides with the expression first derived by Damgaard ||5^. Our 
derivation provides a natural explanation of their result. The two extra fiavors are a 
valence quark and its supersymmetric partner. Remarkably, in the context of chRMT, 
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the spectral density for A not close to zero can still be expressed as the ratio of the same 
two partition functions. It was shown in that (|53D can be reduced to the analytical 
expression for the microscopic spectral density of chRMT which, in the case of zero quark 
masses, is given by [^, ^ 



u 



Psiu) = T^IJnXu) - JNf+liu)JNj-liu)]. 



(54) 



We have thus shown that the microscopic spectral density of chRMT can be derived 
from the partition function of partially quenched chiral perturbation theory. The valence 
quark mass dependence of the chiral condensate can be obtained from (|5^) by replacing 
the sum in (|T3|) by an integral. In terms of the microscopic variables, u = XVHo and 
X = mVEo, this amounts to 

1 , , M 2x 



du- 



■(p(t7^)) 



(55) 



In the thermodynamic limit the spectral density in this equation can be replaced by the 
microscopic spectral density (CT) resulting in 



S(x) 
So 



x[lNf{x)KNf{x) + lNf+l{x)KNf-l{x)]. 



(56) 



By using the duality between flavor and topology EDI the result for the sector of topological 
charge u is obtained by replacing Nj Nf + 

The calculation of the valence quark mass dependence of the chiral condensate from 
the partition function with a noncompact measure would justify both the choice of the 
noncompact measure and the generalization of the arguments from [l^ relating the 
compact and the noncompact measures. In a future publication we will come back to 
this point. In the next section we perform a one-loop calculation of the valence quark mass 
dependence of the chiral condensate. In section 8 it is shown that the zero momentum 
component of this result coincides with the asymptotic expansion of ( p6|) up to this order. 



7. Valence Quark Mass Dependence of the Chiral Condensate in Partially 
Quenched Chiral Perturbation Theory 

From a computational point of view, qChPT and pqChPT are not very different. For 
example, it is much easier to use the quark basis ^ |ET| with mesons qiqj and 

corresponding masses Mfj = Eo(mj + mj)/F'^ (to lowest order), i,j = s,v. To lowest 
order, the bilinear quadratic form in the meson fields is diagonal in the off-diagonal 
mesons fields with quark content q^qj (with i 7^ j). It is therefore trivial to write down 
the propagator for these mesons. However, there are non-trivial interactions between the 



16 



diagonal mesons. Nevertheless it is still possible, as shown by Bernard and Golterman 



45], to invert the inverse propagator analytically. The result for one valence quark mass 



IS 



45 



"^'^^P' p^ + Ml {3 + Nfa){p^ + Mi^{p^ + Mj^){p^ + M^,y > 

where is 1 for the usual mesons and —1 otherwise. The mass of the super-?]' is given 
by: 

. ^f;;/'"f . MiJ±M.^, ,58) 

' l + i\fa/6 l + ]\fa/6 



where y = ml/3M^^ 

The quenched limit is obtained in the limit of infinite sea quark masses at fixed values 
of the valence quark mass and mo. On the other hand, standard ChPT can be recovered 
in two ways. First, in the sector of the mesons with two sea quarks, in the limit where the 
sea quark masses are much less than the mass of the super-?]' particle. In that limit the 
super-?]' can be integrated out of the effective partition and standard chiral perturbation 
theory with Nf flavors is recovered. The propagator of the mesons corresponding to the 
sea quark masses is given by 

^-^""^ = ^TWs ~ ^f^TWs' ^^^^ 

which is just the standard propagator in a quark basis with the singlet channel projected 
out. Second, we consider the limit in which the valence quark masses and the sea quark 
masses are equal and both are much less than the super-?;' mass. Then, the super-?;' mass 
decouples from the theory and, to one-loop order, the contribution from mesons containing 
valence quarks is cancelled by the superpartners. We again recover chiral perturbation 
theory for Nf flavors in both the valence quark sector and the sea quark sector. 

In this work, we are mainly interested in the nonanalytic terms in the quark masses 
contributing to the chiral condensate at one-loop level. We will ignore the corrections 
coming from the O(p^) effective Lagrangian and its unknown effective coupling constants. 
We will consider the limit where the valence quark mass and the sea quark masses are 
much less than the 1]' mass. Therefore, contributions from $o loops proportional to log M^, 
will be ignored as well |^ . 

The thermodynamic limit of the valence quark mass dependence was already calculated 
by Golterman and Leung Here, we consider the generalization of their result to finite 
volume. The volume dependence of the chiral condensate in standard chiral perturbation 



theory has already received a great deal of attention in the literature ^ . We refer to 



these works for technical details. In order to obtain the volume dependence of the quark 
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condensate to one-loop, essentially one has to replace the infinite volume pion propagator 

by 

1 



1 



E 



(60) 



where the sum is over the momenta in the box. In the thermodynamic limit the sum can 
be replaced by an integral resulting in (in four Euclidean dimensions) 

1 A/f2 

where A is the cutoff. 

The generalization of the one-loop result for the valence quark mass dependence of 
the chiral condensate (denoted by a subscript v) to finite volume is then given by 



1 ) NfA(M^ ) — - — •- 

^ ^ '''' 3{1 + Nfy) 



^ ml 



(62) 



Notice that Om^AIM'^) = J2p l/(p^ + M^)^. We study the valence quark mass de- 
pendence of for several special cases. The first limit of interest is the case that both 
the valence quark mass and the sea quark masses are much less than the t]' mass. This 
corresponds to taking the limit y ^ oo of the expression obtained in (|6^) 



So 



1 



In the thermodynamic limit this reduces to 



N'fA{Ml) - A{Ml) + {Ml - Ml) dM^A{Ml) 



(63) 



1 - 



■F'^Nf l^ ("^'' — 2/i ^ 2(m, - 2m„) log — - 



167r' 



(64) 



where fi = A^F^/2So. Notice that the valence quark mass occurs in the operator. Three 
special cases of this result are to be mentioned. For = we recover the one-loop 
result for standard chiral perturbation theory. 

For valence quark masses in the range -C m„ <^ mo the valence quark mass 
dependence of is given by 



,a 2TP4 AT 



(65) 



Below we will show that this result leads to the expression for the slope of the spectral 
density derived in |l56[| . 
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For valence quark masses in the range m„ -C <^ mo the expression 

So 



1 - 



reduces to 
(66) 



The expectation is that on this scale for the valence quark mass the result for should 
reduce to the quenched result. 

Let us consider the quenched limit of (p^). This limit amounts to taking the sea quark 
mesons much heavier than the rj' mass, i.e. the limit y while M^^y = mg/S. In this 
limit (E%) reduces to 



So 
So 



^ 3F^{ 



M2 

vv 



for V GO, 



(67) 
(68) 



where, in the second line, we have only quoted the infrared singular terms. For a valence 
quark masses much less than mo, the term proportional to a can be ignored. 

Let us return to (^). We observe that it corresponds to the quenched result provided 
that mg = 3Mgg/Nf. Indeed, for the low-momentum modes we expect that, according to 
the Witten-Veneziano formula, the topological contribution to the rj' mass is determined 
by the global topological susceptibility. As dictated by the chiral Ward identity, this 
quantity is suppressed by the screening of topological charge due to light sea quarks and 
is given by |^ 



2N. 



f 



(69) 



If we combine this result with the relation between mo and (z/^) given in (^6|) we recover 
the previously mentioned condition. 



8. From pqChPT to chRMT 

Let us study the valence quark mass dependence of in the chiral limit of the sea 
quark masses and valence quark masses in the domain 

Amin < < -pTT • (70) 

^ Aqcd 

As has been argued in section 3, in this domain both the random matrix expressions 
(which follow from the zero momentum component of the partition function (^4]) ) and 
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the pqChPT results ( |62D are valid. For in the range ( [70D the propagator A(M^^) is 
dominated by the zero momentum mode resulting in 

^(^'^) ~ ih- (71) 

The valence quark mass dependence then reduces to 



VMIF"^' 



- So(l-^^). (72) 

In section 6 we derived the valence quark mass dependence of the chiral condensate 
from the zero momentum pqChPT partition function. The result that was obtained by 
integrating the microscopic spectral density coincides with the expression found from 
chRMT |51[|. The asymptotic expansion in powers of 1/x = of T,rmt{x) given 

in (|56| ) with the replacement Nf Nf + |z^| is given by 



^RMT 



M ^ 1 _ + , (^/ + IH + i)(iV/ + IH-|) 



+ 0(-). (73) 



So X 2x2 

In the chiral limit the global topological susceptibility is completely screened. For i> = 
we thus find 

S(x) ~ So (l - , (74) 

which coincides with the result for chiral perturbation theory in this domain. 

In the intermediate region where M^^L ~ 1, pqChPT is still meaningful as a pertur- 
bation theory. The finite volume effects are embedded in the pion propagator. For the 
propagator of a free scalar particle in a finite volume we have EBl 173 



1 A/f^ 

9m.A(M2) = -—^{l + \og—)-g,{M\L), (75) 



where 



167r2 Jo ^0 

and the sum is over a four dimensional lattice of integers. For I/AqcdL"^ <^ M <^ 1/L 
the functions gr can be expanded in powers of M resulting in 

92{M\L) = ]i^ + Y^MM^L^)-| + 0(M^), (77) 
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where (3i = 0.140461 and f32 = —0.020305 are shape coefficients of a 4- dimensional cubic 
box |73|. Let us consider the chiral hmit — * 0. Then, in the intermediate range, the 
asymptotic result (|72D is reproduced. For ML ~ 1 the expansion (|77|) breaks down and 
the functions Qr have to be evaluated numerically. 
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l-N,/mJS„ 
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Fig. 1. The valence quark mass dependence of the chiral condensate S^/Sq as 
a function of the valence quark mass in MeV for two (left) and three (right) 
massless flavors. 



In Fig. 1 we show the ratio S^/Sq as a function of the valence quark mass in MeV for 
two (left) and three (right) massless flavors for a box of size (40 /m)^. The dashed curves 
represent the general expression for the valence quark mass dependence given in (|63|) in 
the limit — >• and A(M^) and 9a,/2A(M^) given by (^) with the exact numerical 
result for the functions Qr- If the numerical results for the Qr are replaced by the leading 
term of the small mass expansion ([77D we obtain the dotted curves (which first lie almost 
exactly underneath the pqChPT and then underneath the chRMT curves). The random 
matrix result for the valence quark masses (|56D for Nf = 2 or Nf = 3 and u = is given 
by the solid curves. We stress that this result can also be derived from the pqChPT 
partition function as shown in section 6. The results in this figure have been calculated 
for A = 770 MeV, F = 93 MeV and Sq/F^ = IQQOMeV. For these values the Thouless 
energy is Ec = F"^ /TjqL? ^ 0.015 Me V^. We observe a region where the random matrix 
results and the pqChPT results are in agreement. The size of this region is consistent 
with the range (P^. For Nf = 2 this region appears larger than for Nf = 3 because the 
leading correction vanishes in the thermodynamic limit for Nf = 2. Because of the large 
box size, the 1/L^ terms are almost invisible in this figure at the point where the nonzero 
momentum modes become important. 
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Fig. 2. The ratio R = {T^rmt — '^q)/{'^v — 5^o) versus the valence quark mass 
m,, for two and three massless flavors. 



In order to study the point where the nonzero momentum modes become important 
we study the ratio 



R 



^RMT 



-En 



(78) 



According to (0) and ([73| ) this ratio is close to 1 in the range (|70D. In Fig. 2 we plot this 
ratio as a function of the valence quark mass for a box size of (40 /m)^ and Nf = 2 
(left) or Nf = 3 (right). At intermediate masses we indeed observe a plateau where 
chRMT and pqChPT are in agreement. Notice that in the center of the overlap region 
([fO|), where m„ scales as the deviation of this ratio from 1 is of order 1/L. Indeed, 

for larger volumes we have observed a more pronounced plateau. The Thouless energies 
can be read off from these pictures. For Nf = 2 and Nf = 3, we find values of around 
0.01-0.02 MeV which is in good agreement with the estimate from = -F^/Sq-L^. 



9. Spectral Density in the Diffusive Domain 

As already has been remarked in section 1, the spectral density is obtained from the 
discontinuity of the valence quark mass dependence of the chiral condensate along the 
imaginary axis. First, notice that the quark mass always occurs in the renormalization 
group invariant combination m{qq) with the sign of (gg) opposite to the sign of m. This 
combination can be rewritten as |m|S. For similar reasons, Sq in the expressions for the 
valence quark mass dependence of the chiral condensate is actually multiplied by \m^\/my 
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with Im-ul = yml. From the elementary identity lime_»o[Y (~3;^+ie) — y (—x^—'je)]/'ix ~ ^ 
we thus conclude that the discontinuity of the first term in (Q) is equal to 2. According 
to ([I^) this results in a contribution to the spectral density of Sq/tt which is in agreement 
with the Banks-Casher formula (|l]). The discontinuities that enter in the calculation of 
the spectral density from the one-loop result for the valence quark mass dependence of 
T,y are given by 

D.cL„„, (^) ^ 2, ,79) 

— - log(m, + |m^|) = log(m^ + x^), (80) 

Tfiy j 

X 

Disc|„^^-^ (log(m^ + \m.o\)) = 2iarctan — , (81) 

TTLg 

DiscL„=„log|m^| =i7r. (82) 

The spectral density can be calculated from the one-loop result for the valence quark mass 
dependence of E^, (|64D by means of {p{\)) /V = Disc|^^^-;^ S„(mt,)/27r. The result is 



V 71 



For rrig ^ the spectral density diverges logarithmically. In the limit A ^ this 
expression reduces to 



(P(A)). 
V 



So 

TT 



1 + 



(iV|-4)Sc 
327r2ArjF4 



lAI 



M) 



which was first obtained by Smilga and Stern ||56 |. 

Our expression for the spectral density can be used to evaluate the connected contri- 
bution to the integrated three-point scalar correlator, 



K 



abc 



d'^xd^yd'^z{0\S''{x)S\y)S''{z)\0), 



^5) 



with S"' = qt"'q and t"" is one of the generators of the flavor SU{Nf) group in the normal- 
ization Tr{t"'t'') = 6"'^/2. This correlator can be expressed in terms of the spectral density 

by m 

^2 - 3A2 



K 



abc 



-mgd' 



abc 



d\ 



[rni + A 



2^3 



(P(A)) 



nis 1 



(86) 



where d"'^'^ are the symmetric structure constants of SU{Nf). By substitution of the 
expression (|83D for the spectral density we then obtain 



V 



-K 



abc 



d' 



AT? - 12 



abc f 



y2 
^0 



NfTUs 1287r2F4' 



(87) 
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which is in complete agreement with the result obtained by Smilga and Stern using stan- 
dard ChPT |5^. Notice that both terms of order 0{m^) and of 0{ml) in (|83|) contribute 



to this result. The existence of contributions of this type was already pointed out in |56 
but partially quenched chiral perturbation theory enabled us to actually calculate them. 

The spectral density in the quenched case can be obtained in exactly the same way by 
calculating the discontinuities of the valence quark mass dependence across the imaginary 
axis. The result is 



(P^(A)) So 

V TT 



. ^0 . \M Spa 



(88) 



A logarithmic divergence of the spectral density in the quenched limit has indeed been 
observed in quenched instanton liquid simulations This result allows us to interpret 
the coefficient of log |A| in the limit A -C of (^) in terms of the screened topological 
charge (p9|). 



10. Conclusions 

We have studied the spectral density of the QCD Dirac operator by means of the 
valence quark mass dependence of the chiral condensate in partially quenched Chiral Per- 
turbation Theory (pqChPT). This is the extension of standard ChPT with an additional 
valence quark. Its mass is required to probe the Dirac spectrum. We have considered two 
different domains of applicability: the ergodic domain with rriy -C F^/L^Sq and the diffu- 
sive domain with F^/L^Sq <^ rriy -C ^qcd- The energy scale F"^ / L'^T.q is the equivalent 
of the Thouless energy in mesoscopic physics. In the ergodic domain the partition func- 
tion is dominated by the zero momentum modes and can be reduced to chiral Random 
Matrix Theory (chRMT). In particular, we have shown this for the microscopic spectral 
density, the valence quark mass dependence of the chiral condensate and the asymptotic 
tail of the two-point spectral correlation function. We have found that a noncompact 
superunitary measure is required to directly reproduce the correct valence quark mass 
dependence. These results strongly suggest that the zero momentum component of the 
pqChPT partition function is equivalent to the chRMT partition function. 

The diffusive regime is the realm of Chiral Perturbation Theory. We have reproduced 
the Smilga-Stern expression for the slope of the spectral density of the QCD Dirac operator 
and have calculated its first order correction in the quark mass. We have shown that these 
results are consistent with standard ChPT. For nonzero sea quark masses the spectral 
density diverges logarithmically with a coefficient that is proportional to the screened 
topological susceptibility. 
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We have studied the transition from chRMT to pqChPT in the domain where both are 
apphcable and have found that both theories give the same valence quark mass dependence 
of the chiral condensate in this domain. The Thouless energy can thus be derived from 
the finite volume scalar propagator. 

If pqChPT provides an accurate description of the QCD partition function we have 
shown that up to the Thouless energy the Dirac spectrum can be described by chRMT. 
Agreement with lattice QCD simulations suggests that this is indeed the case. 



Other scenarios are possible. For example, the proof of the Vafa-Witten theorem 



IS 

not valid for imaginary quark masses. Therefore, breaking of "isospin" symmetry between 
the sea quark sector and the valence quark sector cannot be excluded. Related to this is 
the possibility that the pion decay constant is different for the sea quark channels and the 
valence quark channels. This might restrict the domain of validity of chRMT. However, 
universality arguments suggest that such scenarios are unlikely. Ultimately, these and 
other dynamical questions have to be answered by means of lattice QCD simulations. 
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Appendix A. Calculation of the Microscopic Spectral Density for Nj = 1 

In this appendix we calculate the integrals in the calculation of the microscopic spectral 



density from (?5). We explicitly work out the integrals for Nj = 1 and indicate how they 
can be performed for arbitrary Nf. For Nf = 1 the different factors contributing to the 
integral are given by 

Til o\ (cosgp - COS^i) 

B{cose) = — — , (90) 

(Costco — cos6'2j(cos6'i — costy2j 

^ ^ 2JSo\/K-m,) ^2J^- ^^^^ 
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In the chiral limit rris 0, the determinant can be written as 

gMfcVSocos0i _ gm„yEocoseo _ gm„ VSo cos 0i (92) 

Combining all factors we find 



j=o 8tt^ 





do _ ^iSl 




cos 6*0 — cos 


6*2) (cos ^^1 — COS 6*2) 






2 




^(cos6'o — cos6'i) 



where we have used the anti-symmetry of the measure in Oq and 6i. 
We first perform the integral over 6i 

. f,n |e^^° -e^^^p(cosgi -cosg2) , , . 

'~ J '(cos^^o-cos^Ole^^i -e'^2|2- > 

The integrand can be simplified with the help of the identity ( ^8]) (for s = icj)) resulting 
in 

^ J \l- e-^^oe-i^?i)(e-*^?i _ e^^^^) " ^^^^ 

To evaluate the integral we add opposite infinitesimal imaginary increments to 6*0 and 62 
and take the limit after the integration. For definiteness, let us add a negative imaginary 
increment to ^o- Then we can expand the appropriate terms in a geometric series. We 
observe that 9i only occurs as positive powers of exp(— i^^i), and only the constant term 
contributes. The result for the integral is thus 

h = -2TTe'^^-'^° . (96) 

For the opposite choice of the imaginary increment we obtain the complex conjugate 
expression resulting in the same final answer, 



where x = m^VIlo- Notice that only the real part of Ji contributes to the integral over 
62 and 60 ■ It is simplest to evaluate this integral by first differentiating it with respect to 
X and using once more the identity ( ^81) 

d^j:{x) = y"d^2t^^o(l-e"*^""^^^)(l-e+^^o+^^^)e^^2-^^Oe^(™'^o-™'^^2\ 

^0 

= —[Jii-ix)Ji{-ix) + Jo{-ix)J2{-ix)], 
^0 

= —dx[xJi{—ix)—xJo{—ix)J2{—ix)]. (98) 
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This equation can be integrated trivially with an integration constant that is determined 
by the asymptotic limit for the spectral density. The result for the microscopic spectral 
density is obtained by putting m„ = i\ and taking the microscopic limit according to (§) 
(with u = XVEq), 

Ps{u) = -[J^u) - Mu)J2{u)], (99) 

in complete agreement with the chRMT result for Nf = 1. 

For an arbitrary number of flavors Nf, we systematically use the symmetries of the 
integrand. First, the determinant deto<fc,«<Af/[exp(mfc cos6'i)] appearing in (|45| ) and the 
Berezinian B{cos6i) are both antisymmetric in 6i ^ 6j, for i ^ j and i, j = 0, ■ ■ ■ , Nj. 
Therefore the whole expression can be simplified to: 

P> 7^ f riff M cos 6*0 - COSgjv^+i (^„3ep_^„,g^ ) 

where 



r^,(^^Ar,+i,^^o) = [U dOk |A(e^^^)p cu(0^^+i,0,) ^{-60,6^). (101) 
k=i 

and uj{9,(p) is defined by 

, ^ , , cos 9 — cos d) 

Considering again (92^S(x) and using the fact that |A(e*^'-')p is the sum of 1 plus terms 
exp{iJ2k''^k9k) with finite and J2k''^k = [^], we can show that the only term in 
'TNf{9Nf+i, 9o) contributing to the integral is proportional to exp[iA^/(6'o — ^^at^+i)]. Notice, 
that the term proportional to sin[A^/(6'7v^+i — 6^0)] does not contribute to the integral. 
To obtain the final result for the microscopic spectral density we again use the symme- 
try properties of the integrand with respect to the exchange and translations of the 6k 
variables. 



Appendix B 



In this appendix we calculate the Berezinian corresponding to the diagonalization of 
a superunitary matrix. The calculation is analogous to the case of a super-Hermitean 
matrix []65| , |66| . We can decompose the superunitary matrix U into its eigenvalues A and 
another superunitary matrix V hj U = VAV"^ . Since U is superunitary its eigenvalues 
are just phases, i.e. A = diag(e*^'=). The relevant length element is 



\ds\^ = StrdUdU^ = Str (rfAc/A^) + Str ([vUV,A\ [vUV,A 



(103) 
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Since we are only interested in the A dependence of the Berezinian we can use the 
parametrization V = e^^ so that V'^dV = idH. We then have 



= Y,t{i)dKdK* + ^e(i)e(j)|Ai - Aj\^dHijdHlj, (104) 



where e(z) = 1(— 1) if the Aj corresponds to a fermion (boson). If we write this in the 
form |(isp = Qa^bdzadzl then the Berezinian is given by wSdet(5'), 



n„-|Af-Ap 
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